Metallizing the Mott insulator TiOCl by electron doping: a view from LDA+DMFT 
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Based on recent experiments, we describe the Mott insulating, but undimerized state of TiOCl 
using the local-density approximation combined with multi-orbital dynamical mean field theory 
(LDA+DMFT) for this 3d 1 system. Good agreement with the high estimated value of the superex- 
change is obtained. The possibility of an electron-doped insulator-metal transition in TiOCl is 
investigated in this scheme and a Mott-Hubbard transition with a jump in carrier density inside 
the metallic state is found. Clear, co-incident, discontinuous changes in orbital occupations are 
observed, showing that such a transition would involve strong, multi-orbital correlations. These 
results call for studies on suitably intercalated TiOCl that may induce metallization and, possibly, 
unconventional superconductivity. 
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PACS numbers: 71.27.+a, 71. 15. Mb, 79.60.-i 

Titanium-based oxide systems have recently been 
found to exhibit varied physical properties. Being 
electron-doped analogues (3d 1 ) of the well-known high- 
T c superconductors, they have been studied with a view 
to search for superconductivity. Titanium-oxychloride 
(TiOCl) is of special interest in this context [1]: a 
layered material undergoing spin-Peierls dimerization at 
low T. Unlike other candidate spin-Peierls systems like 
CuGeOs, the superexchange scale J ~ 660 K in TiOCl 
is almost an order-of-magnitude larger. Additionally, the 
fact that the single d electron resides in the t^g orbitals 
might suggest the absence of Jahn- Teller effects, imply- 
ing weakened tendency for doped carriers to localize. It 
has nevertheless proved to be difficult to dope TiOCl us- 
ing traditional techniques. However, structurally similar 
materials such as Li-intercalated ZrNCl are known to 
undergo insulator-metal transition and superconductiv- 
ity [2], leaving the issue open for TiOCl. 

Early observation of spin-gap generation in spin sus- 
ceptibility showing a kink at 94 K and an exponential 
drop below Td — 66 K , accompanied by a doubling in 
the unit cell along b suggested a spin-Peierls dimerized 
state in TiOCl. That this may in fact involve coupled 
spin, orbital and lattice degrees of freedom is suggested 
by large phonon anomalies [3] and by T dependent g- 
factors and ESR linewidths [4]. This should not be too 
surprising: in fact, the magnetism of early TMOs is 
known to be sensitively determined by the above fac- 
tors [5]. Indeed, Seidel et al. [3] and Saha-Dasgupta et 
al. [6] have recently studied the electronic structure of 
TiOCl using the LDA+U method: in contrast to what 
is expected from a cubic crystal field, the low-symmetry 
crystal field in TiOCl splits the t2 g orbitals into a lower 



manages to achieve a spin-Peierls insulating state when 
moderately strong, multi-orbital Coulomb interactions in 
the d shell are switched on is a major question for ab ini- 
tio theoretical approaches. 

In this letter, we aim to address precisely this issue 
in part. Starting with the LDA calculation for the ac- 
tual structure of TiOCl, we explicitly treat the dynami- 
cal effects of strong, multi-orbital electronic correlations 
in the 3d 1 case using multi-orbital dynamical mean field 
theory (DMFT). To solve the DMFT equations, we em- 
ploy the multi-orbital iterated perturbation theory (IPT) 
used successfully earlier [7] . This is in fact necessary: the 
LDA generically gives metallic states owing to neglect of 
strong electron correlations, while LDA+U, which does 
give the correct ground state(s) for insulating systems, 
generically overestimates localization, leading to charge 
gap(s) too large compared to experiment. Moreover, be- 
ing ground state theories, both are incapable of accessing 
the ubiquitous dynamical spectral weight transfers in re- 
sponse to various small perturbations; the importance 
of such aspects is known to be crucial, e.g, in insulator- 
metal transitions. This characteristic underlying the ill- 
understood response of correlated systems can however 
be understood within DMFT and its extensions. 

TiOCl crystallizes in the Pmmn space group at high- 
T, consisting of bilayers of Ti 3+ and 2 ~ parallel to the 
ab plane separated by Cl~ layers. In a perfect octahe- 
dral environment the crystal field would split the atomic 
d levels into three-fold degenerate fa g ) and doubly de- 
generate e g levels. However, the basic TiCliO/^ unit in 
TiOCl is distorted, and the low-symmetry crystal field 
in the real structure splits the t^g degeneracy, with the 
lowest d xy orbital separated from slightly higher lying 



lying d xy and higher-lying d yz ^ zx orbitals. The lone clcc- d yz ^ zx orbitals. The inset of Fig. 1 shows the onc-clcctron 
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tron in TiOCl would then occupy the lowest d xy orbital 
at low-T. However, the higher-lying d yz _ zx orbitals are 
only about 0.2 eV above the lowest one, and, in a situa- 
tion where the LDA bandwidth is about 2.5 eV, strong 
inter-orbital mixing can be expected. How the system 
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density of states (DOS) obtained using density-functional 
calculations in the LDA. As expected, LDA does not give 
an insulating solution. We now show how LDA+DMFT 
reconciles this problem, yielding not only the correct 
ground state, but also a quantitative picture for onc- 
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electron (hole) excitation spectrum. The LDA Hamilto- 
nian is H = Efc, Q/3 ^ c k a a c k0a + Ei,o<7 <w'W> where 
e L<r = ~U{n aS - 1/2) + ^-a(n af7 - 1) with U and 
Jh as defined below. The second term above takes care 
of avoiding the double-counting of interactions already 
treated on the average by LDA. 

We account for strong multi-orbital correlations in 
TiOCl by choosing [6] the intra-orbital U, the local Hund 
coupling Jh and the inter-orbital V ~ U — 2Ju- We 
stress the importance of including U 1 : indeed, neglect- 
ing U' would imply that an electron hopping from site i 
to its neighbor (avoiding the same lowest orbital to es- 
cape the penalty of U) would occupy the slightly higher- 
lying d yz ,zx orbitals. With U' = 0, and a one-electron 
band-width of 2.0 eV, this electron (and the hole it 
leaves behind on the site i) would hop freely, always giv- 
ing a metallic state. The full many-body Hamiltonian 
is H = H + C/Ei,a n i^n iai + Y. iaS 3oo> U^'n iatT ni^. 
For early TMOs, it is sufficient to consider only the 
t 2g bands in the LDA+DMFT scheme [8]. Due to 
tig — e g and tig — 2p hybridization in the real crystal 
structure of TiOCl, we include the TM-e 9 and 0-2p 
parts (of the DOS) with "^g" symmetry in the unper- 
turbed band DOS. Further, in the paramagnetic/ferro- 
orbital situation (see below) applicable to TiOCl in the 
undimcrized state, we have, G a p a(7 > — 5 a pS aa >G aa (uj) 
and Y, a p aa i (w) = 5 a p5 aa iT, a „{(J). 

The DMFT solution in the ti g sector involves (£) re- 
placing the lattice model by a multi-orbital, asymmetric 
Anderson impurity model, along with (ii) a selfconsis- 
tency condition which requires the impurity propagator 
to coincide with the local Green's function of the lat- 
tice, given by G«(w) = ^ /d 3 fc[ ( „ +M)1 _ g o^_ Sn(ij) ]a- 

Using the locality of E Q in d = 00, we have G a (uo) — 
G° a (u> — H a (u>)). Further, U',Jh scatter electrons be- 
tween the d xy and d yZyZX bands, so only the total number, 
E a nf* , is conserved in accord with Luttinger's theorem. 
In the calculations, strict compliance with the Friedel- 
Luttinger sum rule (guaranteeing correct low-energy be- 
havior) and with the first three moments of the spectral 
function (guaranteeing correct high-energy behavior) is 
maintained [9]. 

To solve the multi-orbital single-impurity problem, we 
use the multi-orbital (MO) IPT, suitably generalized for 
arbitrary band-filling and temperatures. This method 
has been successfully used to solve a host of related prob- 
lems [9] , and we refer the reader to the established litera- 
ture for the details. We input the LDA DOS (p° (e)) into 
the DMFT(MO-IPT). In a multi-orbital system, DMFT 
leads to two effects: first, the LDA parameters like crystal 
field splitting are renormalized in non-trivial ways from 
their bare values, and may even change sign [7]. These 
effects result from the multi-orbital Hartree shifts, which 
correspond to effects captured by LDA+U. Secondly, and 
more importantly, DMFT introduces non-trivial effects 



stemming from the dynamical nature of the strong elec- 
tronic correlations. Namely, these processes lead to large 
transfer of spectral weight across large energy scales in 
response to small changes in the bare electronic struc- 
ture, a characteristic lying at the heart of the anomalous 
responses of correlated systems [5]. 
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FIG. 1. LDA+DMFT and LDA (inset) partial densities 
of states (DOS) for TiOCl. Dotted (red) line denotes the 
DOS for the non-degenerate ground-state (xy) orbital, the 
full (blue) line denotes the DOS for the higher (two-fold de- 
generate) lying t2 g orbitals. Higher-lying e g and 0-2p bands 
are not shown. 

Let us now present our results. In Fig. 1, we show the 
orbitally-resolved spectral functions for the d xy (dotted 
line) and the d yZtZX (solid) parts for U — 3.3 eV, with 
Jh = L0 eV and U' — 1.3 eV. As discussed above, 
DMFT results in non-trivial renormalization of the LDA 
crystal field splitting, A t2g = e xy — e a (a = yz, zx) from 
-0.45 eV to 0.12 eV in LDA+DMFT. More important, 
a clear charge gap showing the Mott-Hubbard charac- 
ter of the insulator is selfconsistently generated (the very 
small finite DOS at lj = is a result of finite T « 100 K 
effects). The drastic change in the lineshape attests to 
the explicit consideration of dynamical multi-orbital cor- 
relations. Photoemission spectroscopy should be able 
to confirm this interpretation: if true, a lower Hubbard 
band feature around —1.0 eV with an asymmetric profile 
should be observable. Also XAS should show a shoulder 
feature at low-energies above the Mott gap, and the XAS 
weight should be much higher than in PES. 

An important result directly gleaned from Fig. 1 is 
that Ud xy ■ nd a =7:3, clearly showing that the 
d xy orbital is dominantly populated in the ground state 
(however, with non- negligible inter-orbital fluctuations). 
This is exactly the condition needed to generate one- 
dimensional spin chains along 6, as observed, and is con- 
sistent with the LDA+U result. Using the renormalized 
value of A t2g = 0.12 eV along with the hopping ma- 
trix elements t a p [6] and U.U'.Jh as given above, the 
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effective superexchange constant between nearest neigh- 
bor (S = 1/2) spins on Ti 4+ ions along b can be esti- 
mated. Following [10], we get J ~ 0.06 meV ^ 720 K, 
in nice agreement with the dc susceptibility data. This is 
the dominant contribution to the superexchange, and we 
found that further-neighbor couplings are negligible on 
the scale of J. Given this observation, the spin-Peierls 
instability observed at T$p = 66 K is almost certainly 
caused by coupling to lattice degrees of freedom, and 
not by competing, appreciable next-n.n interactions [6]. 
The derivation of this instability is likely to involve low- 
dimensional effects close to and above T$p and is out of 
scope of any mean-field theory. However, we stress that 
the SP transition arises as an (lattice coupling induced 
spontaneous breaking of translational symmetry along b) 
instability of the (high-T) Mott insulator derived above. 
This can be verified, for example, by studying the evo- 
lution of the dc spin susceptibility upon adding a small 
number of non-magnetic impurities at Ti sites: in anal- 
ogy to Vi- x Cr x C>2 [5], we predict a Curie component 
growing with the impurity concentration, which should 
not be observed if the I-phase is a band-Peierls insulator. 

We now consider the effect of pressure and carrier dop- 
ing: at the outset, we clarify that attempts to dope 
TiOCl have not succeeded to date. On the other 
hand, the fact that related systems like Li x ZrNCl and 
M x HfNCl (M = Li, Na) [2,11] have been shown to be- 
come metallic and superconducting upon electron doping 
suggests that attempts to use similar strategics in TiOCl 
may work [12]. With this caveat, we have studied the is- 
sue of carrier density driven I-M (Mott) transition(s) in 
our LDA+DMFT study. 

Hole doping ( Case 1 ): We have repeated the 
LDA+DMFT calculation for hole-doping (x), corre- 
sponding to band-filling n t2 = 1— x. No instability to the 
metallic state is found up to x = 0.1. We also performed 
the calculations by changing A t2g in small trial steps (be- 
ginning with its value Aj 2g ) with constant n t2g = 1.0, 
as done earlier in another context [7] to mimic pressure 
driven changes, and found no I-M transition. Notice that 
hole doping will predominantly involve changes in the 
lowest, d xy band. The absence of hole-doping driven I- 
M transition is unusual and contrary to naive expecta- 
tions from (d — oo) studies of the doped Hubbard model, 
and its absence under pressure is consistent with high- 
pressure studies which have hitherto found no I-M tran- 
sition up to P — 60 kbar [1]. 

Electron doping ( Case 2): A dramatically different out- 
come is realized with electron {n t2g — 1 + x) doping. A 
nearly first-order I-M transition with rapid change in the 
carrier density around n t2g — 1.9 is clearly seen in the 
LDA+DMFT results. In TiOCl, e~-doping necessarily 
involves doping the d a bands, in contrast to hole doping, 
which affects mainly the d xy band. 

This finding has a natural interpretation in a multi- 
orbital Mott-Hubbard picture. Beginning with Case 2, 
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FIG. 2. t2 g partial DOS for the 3d xz:yz and (inset) 
orbital for different values of the total electron number. 

electrons added to the Mott insulator (in our calculation, 
we change the electron number to n t2g = 1 +x and search 
for the I-M instability starting from the Mott insulator 
found above) prefer to occupy the higher (d a ) orbitals 
to escape the large (U) cost for doubly occupied d xy 
orbitals (U' < U). In this situation, the renormalized 
A t2 = 0.12 (MI) changes upon doping by an amount 
related to U' and to the self-consistently calculated aver- 
age occupations n xy , a of each t2 g orbitals. This in turn 
changes A t2 . More importantly, the dynamical multi- 
orbital correlations (U,U' , Jh) react drastically to this 
small change in A t2g , transferring spectral weight over 
large energy scales from high- to low energies, and driv- 
ing a Mott I-M transition at a critical shown in 
Fig. 2. We clearly observe that only the d a bands show 
metallic behavior; the d xy DOS still represents almost in- 
sulating behavior. This is a direct consequence of the fact 
that it is the d xy orbital which is predominantly occupied 
in the MI. In this limit, the multi-orbital Hartree shift for 
the d a orbitals pulls the d a states lower, populating them 
more (given the constraint), driving more spectral weight 
transfer from high- to low energy (MO-DMFT), sclfcon- 
sistcntly generating a metallic state. The renormalized 
value of A^ g = —0.011 in the metallic (n t 2 g = 1.9) state 
should imply a "melting" of the spin dimerization ob- 
served in the I phase. More importantly, it also implies a 
drastic change in occupation of the various tig orbitals, 
corresponding to the sign-change in A t2 . The situation 
seems to be similar to the one encountered in VO2 [5]. 
In the Ginsburg-Landau picture of the MIT, A t2g plays 
the role of an external field in the orbital sector, which 
drive changes in orbital occupations, leading to the MIT 
via spectral weight transfer (corresponding to the sec- 
ond solution of the DMFT equations) beyond a critical 
x-dependent value. 

Clearly, the case of hole doping is very different: holes 
will predominantly occupy the lowest d xy orbitals. The 
multi-orbital Hartree shift is now ineffective. Since the 
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d xy DOS shows a large (~ 1.0 eV) gap in the /-phase, 
small changes in e xy with x, along with the absence of 
multi-orbital shifts (which produced the lowering of a or- 
bitals in Case 2) are not sufficient to drive an I-M tran- 
sition in Case 1. In reality, given the quasi-lD character 
of the hopping involving d xy orbitals, the doped holes 
would be immediately localized by disorder: this is an 
effect which will operate at least above a certain temper- 
ature scale where the system would really behave one- 
dimensionally. At low-T the above mechanism may be 
operative in reality. 




FIG. 3. Total one-electron spectra for different values of 
the total electron number. The inset shows 3d xz:yz DOS at 
Ef as a function of the the corresponding orbital occupation. 

In Fig. 3, we show the total spectra in the insulat- 
ing and metallic phases. We clearly see the low en- 
ergy pseudogap in the DOS for the bad-metallic case 
(rit 2 < n cr — 1.8) and an abrupt change to a much more 
itinerant (large DOS at low energy) metal for n t2g > n cr . 
Large spectral weight transfer from high energies to the 
Fermi level as a function of electron doping is also clearly 
visible. As remarked earlier, such a dynamical spec- 
tral weight driven I-M transition is out of scope of LDA 
or LDA+U methods, and can only be reliably studied 
within LDA+DMFT. To further explore the non-rigid- 
band scenario across the insulator-metal transition, we 
show the change in p a (Ep) as a function of n a (inset of 
Fig. 3) for our chosen parameter set. Clearly, p a {E F ) 
jumps sharply from 0.17 to 0.31 around n a — 0.56, indi- 
cating an unconventional electron-doped first order tran- 
sition. Our results imply that the metallic state of TiOCl 
is accompanied by an jump in the itinerant carrier den- 
sity at a critical value of the n a orbital occupation close 
to the I-M transition. Based on our calculation, we tenta- 
tively suggest that suitably intercalated (electron doped) 
TiOCl may also exhibit unconventional superconductiv- 
ity upon metallization. Further, given the multi-band 
character of the metallic state found above, multi-band 
(spin singlet for Jh < A t2g or more probably spin triplet 



for Jh > A t2s ) superconductivity may be realized in such 
systems. This issue is out of scope of this paper, however. 

To conclude, we have studied the insulating state in 
TiOCl using the ab-initio LDA+DMFT. In good agree- 
ment with the dc susceptibility measurement, we found 
that the nearest neighbor superexchange, J ~ 720 K 
(it is 660 K in the experiment). A consistent derivation 
of the spin-Peierls instability from the Mott insulator re- 
quires inclusion of additional (non-local?) coupling to the 
lattice, and is planned for the future. In accordance with 
pressure studies, we found no I-M transition occuring un- 
der pressure in TiOCl. Building on this agreement with 
various basic obervations in the I-phase of TiOCl, we 
suggest that electron-doping the Mott insulator (maybe 
achievable by intercalation with suitable species) would 
drive a first-order Mott-Hubbard transition to a corre- 
lated metallic state. Our results show that if this transi- 
tion be first order, it will be accompanied by discontin- 
uous changes in orbital occupations, a feature which is 
also seen in increasing variety of correlated, multi-orbital 
systems undergoing Mott transitions. Our work calls for 
studies on suitably intercalated TiOCl to look for possi- 
ble routes to metallization and unconventional supercon- 
ductivity. 
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